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ABSTRACT 

Observations of jets from young stellar objects reveal the asymmetric outflows 
from some sources. A large set of 2.5D MHD simulations has been carried out for 
axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the 
purpose of assessing the conditions where the outflows or jets are asymmetric rela¬ 
tive to the equatorial plane. We consider initial magnetic fields that are symmetric 
about the equatorial plane and consist of a radially distributed field threading the 
disc (disc-field) and a stellar dipole field. (1). For pure disc-fields the symmetry or 
asymmetry of the outflows is affected by the midplane plasma /3 of the disc (where 
P is the ratio of the plasma pressure to the magnetic pressure). For the low density 
discs with small plasma P values, outflows are observed to be symmetric to within 
10% over timescales of hundreds of inner disc orbits. For the denser higher [3 discs, 
the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic 
field motion in the two hemispheres becomes different. This asymmetry leads to 
asymmetric episodic outflows. ( 2 .) Accreting stars with a stellar dipole field and no 
disc-field exhibit episodic, two component outflows - a magnetospheric wind and an 
inner disc wind from somewhat larger radial distances. Both are characterized by 
similar velocity profiles but the magnetospheric wind has densities > 10 times that 
of the disc wind. (3.) Adding a disc-field parallel to the stellar dipole field acts to en¬ 
hance the magnetospheric winds but suppress the disc wind. (4.) In contrast, adding 
a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress 
the magnetospheric and disc winds. Our simulations reproduce some key features of 
observations of asymmetric outflows of T Tauri stars. 

Key words: accretion, accretion discs - MHD - black hole physics, magnetic fields, 
jets, stars: winds, outflows 


Astrophysical jets are observed from many disc accreting 
objects ranging from Young Stellar Objects (YSOs) to sys¬ 
tems with white dwarfs, neutron stars and black holes 
(BHs) (e.g., Livio 1997). The jets transport angular mo¬ 
mentum away from the disc and thereby facilitate accretion 
to the central object. A large-scale magnetic field thread¬ 
ing the disc is thought to have a key role in launching and 
collimating the jets (e.g., Lovelace et al. 2014). 


1 INTRODUCTION 


There is clear evidence mainly from Hubble Space Tele¬ 


scope (HST) observations that YSOs often show an asym¬ 
metry between the jet and the counter jet. In HH30 the 
jet has more structure than the counter jet in the form of 
knots, suggesting different mass densities or different accel¬ 
eration mechanisms (Mundt et al. 1990; Ray et al. 1996). 
The mass flux of the jet is approximately constant along its 
entire length whereas the counterjet’s mass flux varies by a 
factor of 2 along its length (Bacciotti et al. 1999). RW Aur 
A exhibits a jet with velocity approximately twice that of its 
counter jet and a mass flux 2 — 3 times as large (Woitas et al. 
2002, Hartigan & Hillenbrand 2009, Melnikov et al. 2009). 
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Likewise LkHa 233 exhibits a jet with a knotted structure, 
whereas its counter jet has no knots but has similar velocity 
(Perrin & Graham 2007). FS Tan B has been observed to 
have a jet with a higher velocity but a lower density than 
its counter jet (Liu et ah 2012). Rotation measurements of 
T Tauri jets hnd evidence of jets counter-rotating with re¬ 
spect to the accretion disc or of the opposite rotation of 
jets in the northern and southern hemispheres (Sauty et al. 
2012). Studies have shown that the jet/counter jet structure 
are different, with radial velocity differences of ~ 10 — 25 
km/s (Coffey et al. 2004; Ray et al. 2007 for a review). 
One of the few exceptions is HH 212 which exhibits nearly 
symmetric jet/counterjet (Zinnecker 1998). Ellerbroek et al. 
(2013) found HH 1042 to have jets with similar mean veloc¬ 
ities but the knotted structure and velocity spread suggests 
an asymmetric launching mechanism. Whelan et al. (2014) 
have studied the jet/counter jets mass outflow rates and 
stellar accretion rates in ESO-Ha Par-Lup 3-4 to constrain 
lauching mechanisms. Though most observations show some 
sort of asymmetry in the outflows, this may be a result of 
asymmetries in the ambient medium and not the launching 
mechanism (Carratti o Garatti 2013, White 2014). 

Relatively little theoretical and computational work 
has been directed at understanding the asymmetry of the 
YSO jets. A theoretical study by Wang et al. (1992) 
considered disc magnetic fields composed of dipole and 
quadrupole components and found that the magnetically 
driven jets could in some cases be highly asymmetric. 
Ghagelishvili et al. (1996) investigated how a quadrupo- 
lar disc dynamo interacting with a stellar dipole could give 
rise to one sided jets. Lovelace et al. (2010) carried out 
axisymmetric MHD simulations of disc accretion to stars 
with magnetic field consisting of different combinations of 
aligned dipole and quadrupole field components and with¬ 
out the assumption of symmetry about the equatorial plane. 
The simulations showed that for conditions where there is a 
significant quadrupole component in addition to the dipole 
component, a dominantly one-sided conical wind tends to 
form on the side of the equatorial plane with the larger value 
of the intrinsic axial magnetic field. Eor cases where the 
quadrupole component is absent or very small, dominantly 
one-sided outflows also form, but the direction of the flow 
flip-flops between upward and downward on a time-scale 
of ~ 30 d for a protostar, von Rekowski (2003) and von 
Rekowski & Brandenburg (2006) investigated outflows gen¬ 
erated by time dependent stellar and disc dynamos. Mat- 
sakos et al. (2012) investigated how field alignment between 
the disc and star and pressure distribution in the ambi¬ 
ent cloud caused velocity asymmetries in protestellar jets. 
Eendt and Sheikhnezami (2013) investigated how thermal 
and density asymmetries in the disc can lead to asymme¬ 
tries in the outflow. Velazquez (2014) performed 3D hydro- 
dynamical simulations to study the case of asymmetric jet 
launching from protoplanetary discs that are preferentially 
fed matter from one side of the disc by a binary companion. 
Stute et al. (2014) showed that relaxing axisymmetry has 
little effect on MHD driven winds. 

Glassical T Tauri stars (GTTS) commonly show evi¬ 
dence of a stellar dipole and higher order multipoles with 
held strengths ~ 1 — 2 kG at the stellar surface (Jardine, 
Gollier Gameron V Donati 2002; Donati et al. 2007, 2008). 


Eurthermore, it is known from MHD simulations that any 
large-scale magnetic held threading the star’s accretion disc 
tends to be advected inward due partly to viscous trans¬ 
port in the disc and partly due to angular momentum loss 
to outhows (e.g., Dyda et al. 2012). A large-scale magnetic 
held is thought to be responsible for driving and collimat¬ 
ing the jets from young stellar objects. The origin of the 
required magnetic helds is not established, though a strong 
possibility is that dynamo processes in the rotating, convec¬ 
tive stars. Eor example, the Zeeman signature in GQ Lup 
(Donati et al. 2012) and DN Tan (Donati et al. 2013) sug¬ 
gests that the magnetic helds of T Tauri stars are generated 
by non-stationary dynamos. Observations of disc-accreting 
stars show strong correlation between the magnitude of the 
star’s dipole magnetic held and the rotation period of the 
star (Vidotto et al, 2014). It is commonly thought that the 
observed magnetic helds of young stars arise from dynamo 
processes inside the rotating convective stars. However, it 
is possible that dynamo processes operate in both the ac¬ 
cretion disc and in the star as considered by von Rekowski 
and Brandenburg (2006). One then has the possibility that 
the magnetic held generated in the disc is accreted onto 
the star and thereby modihes the stellar held changing the 
star’s dipole, quadrupole, and octupole moments for ex¬ 
ample (Dyda et al. 2015). Time dependent dynamo pro¬ 
cesses could also be responsible for driving periodic or quasi- 
periodic outhows (Stepanovs, Eendt & Sheikhnezami 2014). 
It is therefore important to understand the interaction of 
the dipole held of a rotating star with the disc held. In 
particular, we are interested in how the symmetry of the 
stellar plus disc held conhguration affects the symmetry or 
asymmetry of the out how or jets. 

We consider four classes of magnetic held conhgura- 
tions - a pure disc held (DISG), a pure stellar dipole held 
(DIP), a disc held with anti-aligned stellar dipole held 
(ADIP) and a disc held with parallel aligned stellar dipole 
held (PDIP). 

In DISG runs we consider an initial reference disc-field 
the same as that studied by Zanni et al. (2007) in axisym¬ 
metric PLASH code simulations which assumed symmetry 
about the disc plane. These simulations exhibited a sta¬ 
ble, symmetric inner disc wind. Later work by Eendt and 
Sheikhnezami (2013) considered the same disc held but 
without assuming symmetry about the equatorial plane. 
They found that this held gave rise to symmetric out¬ 
hows, but that asymmetric outhows occurred if the disc 
density and temperature were asymmetric about the equa¬ 
torial plane. 

In addition to the mentioned disc-held, we consider 
accretion onto magnetized rotating stars with an aligned 
dipole magnetic held with magnetic moment /x. Previous 
work by Lovelace et al. (2010) found that a pure dipole 
in the “propeller regime” exhibits asymmetric outhows. In 
the propeller regime, the radius of the star’s magnetosphere 
Rm (where the pressure of the disc matter is of the order 
of the magnetic energy density of the star’s dipole held) is 
larger than the corotation radius Rcr = (where 

the rotation rate of the star equals the Keplerian rotation 
rate of the disc) (Illarionov & Sunyaev 1975). In the pro¬ 
peller regime a large fraction of the accreting disc matter 
may be expelled from the system (Lovelace, Romanova V 
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Bisnovatyi-Kogan 1999; Romanova et al. 2004, 2005; Ustyu- 
gova et al. 2006; Lii et al. 2014). A qualitative picture of the 
accretion is that matter accumulates at the magnetospheric 
boundary and compresses the magnetosphere inward until 
Rm < Rcr and funnel flow accretion occurs. The funnel 
flow spontaneously breaks the system symmetry about the 
equatorial plane and strong outflows are generated in the 
opposite hemisphere (Lii et al. 2014). As the matter reser¬ 
voir at the magnetospheric boundary is depleted, the mag¬ 
netosphere expands outwards and the system can “forget” 
which side had a funnel flow. This process repeats on time 
scales of roughly 40 days for CTTS, with the side deter¬ 
mining the funnel flow being randomly determined, leading 
to episodic, asymmetric outflows (Lii et al. 2014). The DIP 
simulations are used to compare our simulations with stellar 
dipoles to these previous works. 

Moreover, we consider the two cases where the stellar 
dipole field at the magneotspheric radius is parallel (PDIP) 
or anti-parallel (ADIP) to the disc-field. This offers two 
types of behaviors which can break the symmetry about 
the equatorial plane. In the anti-parallel case (ADIP), the 
stellar dipole field may reconnect with the disc field as it 
is advected inward. Analogous reconnection occurs between 
the solar wind magnetic field and the magnetic field of the 
earth when the two fields are anti-parallel (e.g., Hargreaves 
1992). This reconnection occurs at the magnetopause (in¬ 
side the Earth’s bow shock) where the subsonic solar wind 
runs inot the Earth’s magnetosphere. In the parallel case 
(PDIP), the disc field causes additional pressure at the mag¬ 
netospheric boundary, reducing the magnetospheric radius 
and altering the field line geometry outside the magneto¬ 
sphere. 

In Sec. 2 we discuss the numerical methods used and 
the simulation parameters. Section 3 discusses the four main 
cases simulated. Section 3.1 describes results for the pure 
disc-field (DISC) for different disc densities or plasma /3 
values. Section 3.2 describes results for a pure stellar dipole 
field (DIP), a disc-field and a parallel stellar dipole (PDIP) 
and a disc-field and an anti-parallel stellar dipole (ADIP). 
Section 4 gives the main conclusions from this work and its 
implications for observations. 


2 THEORY 

2.1 Basic Equations 

The plasma flows are assumed to be described by the equa¬ 
tions of non-relativistic magnetohydrodynamics (MHD). In 
a non-rotating reference frame the equations are 


|+V.(,v) = 0 . 

(la) 


(lb) 

<9R 

^+cVxE = 0, 

(Ic) 

^^+V-(pv5) = Q. 

(Id) 
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Here, p is the mass density, S is the specific entropy, v is the 
flow velocity, B is the magnetic field, T is the momentum 
flux density tensor, E = — vxB/c+? 7 tVxB/c is the electric 
field, Q is the rate of change of entropy per unit volume due 
to viscous and Ohmic heating in the disc, c is the speed of 
light and pt is the turbulent magnetic diffusivity. We assume 
that the heating is offset by radiative cooling so that Q = 0 . 
Also, g = — (GM/r^) r is the gravitational acceleration due 
to the central mass M. We model the plasma as a non- 
relativistic ideal gas with equation of state 


= In 

\P^ 


( 2 ) 


where p is the pressure and 7 = 5/3. 

In most of this paper we use cylindrical coordinates 
(R, 0, Z) as these are the coordinates used by our code. 
However, for some purposes spherical coordinates are ad¬ 
vantageous, and they are denoted (r, ^, 0 ). 

Both the viscosity and the magnetic diffusivity of the 
disc plasma are thought to be due to turbulent fluctuations 
of the velocity and magnetic field. Outside of the disc, the 
plasma is considered ideal with negligible viscosity and dif¬ 
fusivity. The turbulent coefficients are parameterized using 
the a-model of Shakura and Sunyaev (1973). The turbulent 
kinematic viscosity is 


JAt 




( 3 ) 


where Cg is the midplane sound speed, Qk is the Keple- 
rian angular velocity at the given radii and ai, ^ I is a 
dimensionless constant. Similarly, the turbulent magnetic 
diffusivity is 


Pt = 


Qk 


( 4 ) 


where an is another dimensionless constant. The ratio. 


r = 


( 5 ) 


is the magnetic Prandtl number of the turbulence in the disc 
which is expected to be of order unity (Bisnovatyi-Kogan & 
Ruzmaikin 1976). Shearing box simulations of MRI driven 
MHD turbulence in discs indicate that V ^ 1 (Guan & 
Gammie 2009, Lesur & Longaretti 2009) and this condition 
has been used in previous disc wind studies (see for example 
Tzeferacos et al. 2013). All our simulations use = ar^ = 
O.I. 

The momentum flux density tensor is given by 


'T - A ^ 

lik — P^ik T pViVk ~\~ I „ Oik . 

' Stt 47r 


+ Tik 


( 6 ) 


where Tik is the viscous stress contribution from the tur¬ 
bulent fluctuations of the velocity and magnetic field. As 
mentioned we assume that these can be represented in the 
same way as the collisional viscosity by substitution of the 
turbulent viscosity. Moreover, under most conditions found 
in our simulations \vz\ ^ ^ Vf) within the disc. There¬ 

fore we assume that the viscous stress is dominated mainly 
by terms involving the angular velocity gradient and the ra¬ 
dial accretion velocity and choose to neglect terms involving 
vz- The leading order contribution to the momentum flux 
density from turbulence is therefore 
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Parameters 

Symbol 

Value 

mass 

Mo 

1.59 X 10^3 g 

length 

Ro 

1.50 X lO^^cm 

magnetic field 

Bo 

8.04 X 10-2 G 

time 

Po 

3.53 X 10-2 y 

velocity 

VO 

8.42 X lO^cm/s 

density 

po 

7.24 X 10-15 g/cm3 

accretion rate 

Mo 

2.72 X lO-i Mo/yr 

temperature 

To 

4.26 X lOl K 

dipole strength 

p 

2.69 X 1035 Q emS 


Table 1. Mass, length, and magnetic field scales of interest and 
the corresponding scales of other derived quantities for a CTTS. 
One can obtain the physical values from the simulation values by 
multiplying by the corresponding dimensional quantity above. 


TR(f) — 

rzR — 

T(j)^ = 


dvR 

-2u,p- , 


TZ4) — 
trr — 


„ dvR 


( 7 ) 


where O = is the plasma angular velocity. 

The transition from the viscous-diffusive disc to the 
ideal plasma corona is handled by multiplying the viscosity 
and diffusivity by a dimensionless factor ^(p) which varies 
smoothly from ^ \ iox p ^ pd — 0.75p(i?, Z = 0) to 

^ = 0 for p ^ 0.25p(i as described in Appendix B of Lii, 
Romanova, & Lovelace (2012). 

The MHD equations are written in dimensionless form 
so that the simulation results can be applied to different 
types of stars. We assume that the central object is a CTTS 
with mass M* = O.SM©, a radius R* = and a mag¬ 
netic field magnitude R* = 3 x 10^ G which is typical for 
the magnitude of the stellar dipole. We define dimension¬ 
ful quantities with a 0 subscript, to denote typical values 
of physical parameters at a reference radius Ro- The mass 
scale is set by the stellar mass so Mo = M*. The reference 
length, Ro = O.IAU, is taken to be the scale of a typical 
inner disc radius. Assuming a stellar dipole field, the mag¬ 
netic field strength Bq — jjR qY. The mass, length 
and magnetic field scales allow us to define all other dimen¬ 
sionful quantities. 

The reference value for the velocity is the Keplerian 
velocity at the radius Ro, '^’o = (GMo/Ro)^^^- The ref¬ 
erence time-scale is the period of rotation at Ro, Po = 
^ttRq/vq. From the MHD equations, we get the relation 
PqVq — Ho/47r, where po is a reference density. This allows 
us to define the reference density po of the disc. The refer- 
ence mass accretion rate is Mo = A'kpqVqRq. The reference 
temperature Tq — Vq/2x mnlkB where ttih is the atomic 
mass of hydrogen and ks is the Boltzmann constant is ob¬ 
tained by taking the ratio of the reference pressure and 
reference density. One can obtain dimensionful quantities 
from simulation results by multiplying by the appropriate 
dimensionful quantity in Table 

Results obtained in dimensionless form can be applied 


to objects with widely different sizes and masses. However, 
the present work focuses on CTTS with the typical val¬ 
ues shown in Table A more detailed description of our 
code and numerical methods can be found in Koldoba et 
al. (2015). 


2.2 Initial Conditions 


2.2.1 Magnetic Field 


The initial magnetic fields are taken to be linear combina¬ 
tions of the disc-field (Zanni et al. 2007) defined by the flux 
function 


'i>d = iBoRt ( 


3/4 


, 5/4 


RoJ {m? + ’ 


( 8 ) 


and a fixed dipole field centered in the star which has the 
flux function 




( 9 ) 


where p = const is the star’s dipole moment. Our base 
case is taken to be p = 30, but we also consider cases with 
p = 10, 20. The dimensionless parameter m in Eqn. 8 de¬ 
termines the bending of the disc-field lines in the poloidal 
plane with the limit m ^ oo corresponding to purely ver¬ 
tical field lines. We assume m = 1. The poloidal magnetic 
field components are given by 


^ RdR' 


Br = - 


1 dT 

R^ ■ 


( 10 ) 


The relative strengths of the disc and dipole field compo¬ 
nents are characterized by the radius Rr where the stellar 
dipole and disc-field have equal values in the midplane of 
the disc at t = 0. We take Rr = 2.8,4.1, & 5. The magne- 
tospheric radius Rm is determined by plotting the angular 
velocity of the matter in the equatorial plane as a func¬ 
tion of R. The change from the constant stellar rotation to 
Keplerian occurs at Rm- As mentioned in §1, the propeller 
regime occurs for Rm > Rev 

Our simulations cover a wide range of dipole and disc- 
field combinations. In particular, we consider a pure disc- 
field (DISC), a pure stellar dipole (DIP), a dipole and par¬ 
allel disc-field (PDIP), and a dipole with an anti-parallel 
disc-field (ADIP). Figure shows the initial disc density 
and poloidal field lines for the disc-field DISC. Figure 
shows the same quantities for cases including an aligned 
stellar dipole magnetic field (DIP, PDIP & ADIP). 


2.2.2 Matter Distribution 


Initially, matter of the disc and corona are assumed to be in 
mechanical equilibrium (Romanova et al. 2002). The initial 
density distribution is taken to be barotropic with 


pip) 


p/Tdisc, p>Pb and R^Rd, 
p/Tcoi, p < Pb or R^ Rd , 


( 11 ) 


where pb is the level surface of pressure that separates the 
cold matter of the disc from the hot matter of the corona 
and Rd is the inital inner disc radius. The initial inner radius 
of the disc is Rd = 5Ro. At this surface the density has 
an initial step discontinuity from value p/Tdisc to p/Tcov 
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Figure 1. The initial density p and poloidal field lines for radially distributed disc-field in the cases of Bernoulli parameters Left: 
Ho = 3 X 10 Center: Bq = 5 x 10 and Right: Bq = 1 X 10 Increasing the Bernoulli parameter increases the midplane disc 
density and the disc scale height. 



The initial dimensionless temperature in the disc is Tdisc = 
(p/p)disc = 5 X 10“^, and the initial temperature in the 
corona is Tcor = {p/p)cor = 0.5. 

Because the density distribution is barotropic, the ini¬ 
tial angular velocity is a constant on coaxial cylindrical sur¬ 
faces about the z—axis. Consequently, the pressure can be 
determined from the Bernoulli equation 

F(p) + $ + $e = i3o , (12) 

where Bo is Bernoulli’s constant,$ = —GM/\/B? + is 
the gravitational potential with GM = 1 in the code, = 
rdr Q?(r) is the centrifugal potential, and 

Fi ) = f ^disc ln(p/p6) p>ph and Rd , 

^ Tcor \n{p/pb) p <Pb or R^ Rd . 

Varying Bo allows one to vary the initial midplane disc den¬ 
sity. We consider three values of the Bernoulli parameter 
i3o = 3 X 10“^, 5 X 10“^ and 10“^, which we refer to as the 
diffuse, intermediate and dense disc cases respectively. The 
initial half-thickness of the disc h is taken to he h/R = 0.100 
at the inner disc radius Rd — 'o. 

Figure shows the radial dependence of the midplane 
disc density for the different Bernoulli parameters. Fig¬ 
ure shows the initial radial dependence of the midplane 


13 = S7vpIB‘^ (the ratio of the gas pressure to the magnetic 
pressure) for the different Bernoulli parameters. 

2.2.3 Angular Velocity 

The initial angular velocity of the disc is slightly sub- 
Keplerian, 

Q = {1 - 0.003)Qk{R) R>Rd, (14) 

Inside of Rd, the matter rotates rigidly with angular velocity 

Q = {l- 0.003)QK{Rd) R^Rd- (15) 

The corotation radius Rev is the radius where the angular 
velocity of the disc equals that of the star; that is, Rcr = 
{GM/Q,^)^^^. In this study we have chosen Rcr = 1.5. For 
the considered stellar dipole fields and accretion rates, the 
systems are in the strong propeller regime where Rm is 
significantly larger than Rev 

2.3 Boundary Conditions 

Our simulation region has three boundaries: the axis, the 
surface of the star and the external boundaries. For each 
dynamical variable we impose a boundary condition consis¬ 
tent with our physical assumptions. 
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We assume axisymmetry about the axis. On the star 
and the external boundaries we want to allow fluxes and 
impose free boundary conditions dTjdn — 0 where T is 
a dynamical variable and n is the vector normal to the 
boundary. 

At the external boundary along the edge of the disc 
we allow new matter to flow into the simulation region. We 
impose the condition that the matter must be accreting 
Ur < 0. In the coronal region, we prescribe outflow condi¬ 
tions and allow matter, entropy and magnetic flux to exit 
the simulation region. 

The star is taken to be cylindrical in shape with radius 
= 1 and height — 2. Free boundary conditions are 
imposed on the variables p, p and B. The velocity vectors 
are required to align with the magnetic held lines on the star 
via V = (B'v)B/|B|. In addition, we require that matter on 
the stellar boundary how into the star. We treat the corner 


of the star by averaging over the nearest neighbour cells in 
the R and Z directions. 

In addition to the boundary conditions imposed for the 
well-posedness of our problem we impose additional condi¬ 
tions on variables at the boundaries to eliminate numerical 
artifacts of the simulations. 

Our simulation uses a grid of 154 x 225 cells. The star 
has a radius of 1 in units of the simulation and is cylindri¬ 
cal in shape. It extends 10 grid cells above and below the 
equatorial plane. In the R-direction, the hrst 30 grid cells 
have length dR = 0.05. Later cell lengths are given recur¬ 
sively by dRi^i = 1.025dRi. Similarly, in the Z-direction 
the hrst 30 grid cells above and below the equatorial plane 
have length dZ — 0.05. Later cell lengths are given recur¬ 
sively by dZj-^i = 1.025dZj. Figurej^shows a sparse version 
of our grid with only every 7th gridline in the R—direction 
and every 9th gridline in the Z—direction are shown. 



R 


Figure 3. Initial midplane density prohle p(R, 0) as a function 
of radius for different values of the Bernoulli parameter Bq = 
3 X 10-4,5 X 10-4 and IQ-^. 



Figure 4. Initial radial dependences of the midplane (3 = 
Sttp/B^ for the different Bernoulli parameters. 
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3 RESULTS 

We are interested in the mass, angular momentum, and en¬ 
ergy outflow from the disc. We have also calculated these 
fluxes at larger values of |Z| and found this did not signif¬ 
icantly change the values. We argue it is best to avoid the 
last few simulation region cells to avoid boundary effects 
while being as far away from the star as possible to allow 
the matter to be accelerated by the magnetic field. We note 
that mass fluxes are more symmetric when measured near 
the surface of the disc but this is observationally less inter¬ 
esting than the flux at large Z. In addition, for the mass 
flux we require that Vz > 0.1 which corresponds to an ob¬ 
servationally interesting value of > 200km/s. A summary 
of our runs, their key parameters and results are found in 
Table [a 


3.1 Pure Disc-Field DISC 

In this series of runs we varied the Bernoulli parameter 
Bo which determines the initial matter configuration (Eqn. 
12). We take Bo = 3 x 10“^, 5 x 10“^ and 10“^. Increasing 
the Bernoulli parameter increases the midplane density, and 
therefore the total mass of the disc. The initial midplane 
density profiles are shown in Fig.[a For simplicity, we refer 
to them as the diffuse or low-/3 (Bo = 3x 10“^), intermediate 
(Bo = 5 X 10“^), and dense or high-^d (Bo = 10“^) discs, 
respectively. 

For the same disc-fields, the three disc configurations 
exhibit qualitatively different behaviors of their outflows. 
The diffuse low-/3 disc exhibits symmetric (=LZ), long last¬ 
ing outflows as observed by Fendt and Sheikhnezami (2013). 
The intermediate disc exhibits oscillatory behaviour, with 
the jet flip-flopping several times before settling at late time 
to steady symmetric outflows. The dense high-/3 disc ex¬ 
hibits outflows which are persistently asymmetric. This is 


Type 


Bo 


Outflow Summary 


0 

1 

3 

Long lasting, symmetric 

DISC 

0 

1 

5 

Quasi-periodic 


0 

1 

10 

Persistently asymmetric 


10 

0 

3 

Strong disc/mag. wind 

DIP 

20 

0 

3 

Strong disc/mag. wind 


30 

0 

3 

Strong disc/mag. wind 


10 

1 

3 

No magnetospheric wind 

PDIP 

20 

1 

3 

Weak disc/strong mag. wind 


30 

1 

3 

Weak disc/strong mag. wind 


10 

-1 

3 

Suppressed 

ADIP 

20 

-1 

3 

Suppressed 


30 

-1 

3 

Suppressed 


Table 2. Summary of key parameters in our runs and a quali¬ 
tative summary of the outflows observed. 

the result of the excess accumulation of matter in one or 
the other of the two hemispheres near the magnetosphere 
which then feeds more matter into one or the other out¬ 
flows. Murphy, Ferreira and Zanni (2010) have shown that 
the important criteria for launching outflows is that /3 ~ 1 
at the disc surface. Our results show that the character of 
these outflows is affected by the midplane /3. 

d.i.i Diffuse Low-[3 Disc 

Figure shows the mass flux M through the upper hemi¬ 
sphere (solid red line) and lower hemisphere (dashed blue 
line) of the low [3 disc as a function of time. This configu¬ 
ration exhibits outflows with mass fluxes of order 0.1% of 
the stellar accretion rate. For ease in comparing the out- 



R 


Figure 5. Sparse version of grid used in the simulations, showing 
every 7th gridline in the direction and every 9th gridline in 
the Z—direction. 
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Figure 6. Diffuse low—^ disc: The lower plot shows the mass 
flux M through the upper (solid red line) and lower (dashed blue 
line) boundaries. The middle plot shows the mass flux difference, 
normalized to average total mass ffux AM = (Mt — Mi))/{Mt + 
Mf) as a function of time. The upper plot shows the average 
velocity of outffows, weighted by mass along the upper (red) 
and lower (blue) boundaries. Both the mass ffux and velocity of 
outffows is symmetric in this case. 


flows, we plot the average mass flux asymmetry AM = 
{Mt — Mh)/{Mt + Mb) as a function of time. We notice that 
t < 400 is characterized by sporadic asymmetries of order 
10% of the outflow rate. For t > 400 the asymmetry is still 
roughly 10% but is persistently in the lower half plane. The 
upper contour plot shows the average velocity, weighted by 
mass, of outflows along the upper (red) and lower (blue) 
boundaries. 

This serves as a calibration run in the respect that it 
reproduces the results of Fendt and Sheikhnezami (2013) 
which found persistent, symmetric outflows for t < 500 
in our simulation units. Their outflow was measured near 
the surface of the disc, at a distance oi z — 3H. We note 
that our outflows are more symmetric when measured at 
smaller heights, notably the small asymmetry near t ~ 100 
vanishes. However, as we are ultimately interested in out¬ 
flows which exceed the local escape velocity, we choose to 
measure outflows closer to the outer boundary to allow the 
outflows to be magnetically accelerated. We also note that 
both outflows are launched with comparable velocities. 


3.1.2 Intermediate Disc 

Figure shows the mass flux M through the upper hemi¬ 
sphere (solid red line) and lower hemisphere (dashed blue 
line) as a function of time. After an initial period of sym¬ 
metric outflows, we observe for times 100 < t < 400 quasi- 
periodic, oscillating, outflows. For times t > 400 the out¬ 
flows are again symmetric with asymmetries on the order 
of a few percent. 

The oscillations in the outflows are caused by a peri¬ 
odic motion of the field lines. In one hemisphere the field 


lines become more vertical, reducing the outflow, whereas 
in the other they become more inclined relative to the 
Z—axis. This increases the outflow, as matter is more read- 
illy launched along inclined field lines (e.g., Lii et al. 2012). 
As shown in Fig. the structure then reverses, the more 
vertical field lines becoming more horizontal and vice versa. 
This leads to quasi-periodic oscillations in the mass outflows 
with a period T ~ 80 and an amplitude of roughly 50%. At 
the peak of the oscillation, there is therefore a factor of 
~ 3 between the upper and lower mass flux. For a T Tauri 
star with our choice of reference units this corresponds to 
variations on a time scale of roughly 2.8 years. The oscil¬ 
lations are caused by Alfven waves which propagate with 
a speed va ~ ~ The time-scale r for 

passage of the wave across the disc a distance ~ 2h is 
T ~ 2\f^R/vK ~ . We find a midplane fd ^ 4D 

and take ~ 3, the radius from which strong outflows are 
launched and find r ~ 60, which agrees approximately with 
the observed period T ~ 80. 

We note that though the peak mass flux through each 
hemisphere is comparable during this oscillatory phase, the 
velocity of the outflows is not. In particular, the bottom 
jets velocity vz ~ —0.1 during mass flux minima but rises 
to Vz ~ —0.8 and —1.5 during its two peaks. The top jet 
on the otherhand has vz ~ 0.2 during outflow minima and 
peaks at vz — 0.3 and 0.4. A mass flux maxima of the 
bottom jet therefore could explain RW Aur A observations 
for instance, where the jet has a mass flux of 2-3 times the 
counter jet and a velocity roughly twice as great. 

At late times, t > 600, a persistent asymmetry of the 
outflows is observed. This occurs because the original mat¬ 
ter in the main part of the disc has been depleted and new 
matter is being added along the right boundary. This in¬ 
duces an asymmetry in the outer part of the disc which 
in turn induce disc oscillations that lead to a tilted disc. 
This asymmetry appears to be due to our outer boundary 
condition. 


3.1.3 Dense High-13 Disc 

A more massive disc results in a higher density of matter 
close to the star. At late times the disc half thickness h/R 
is roughly 5 times greater than for the diffuse disc in the 
very inner disc R < 5. The matter available for outflows is 
hence greater but is also more strongly influenced by inner 
disc oscillations. Radially propagating disc waves matter 
preferentialy accumulates in the upper or lower hemisphere 
leading to stronger outflows. However, unlike in the dipole 
case where the magnetosphere re-expands after a strong 
outflow, reestablishing a symmetric disc, these asymmetries 
persist until strong vertical disc oscillations develop or the 
asymmetric matter reservoir is depleted. 

As in the case of the intermediate disc, there are strong 
oscillations of the inner disc held. These lead to spontaneous 
one sided outflows as the held lines become more inclined 
to the Z—axis (Fig. |^. The outflows turn off as the held 
becomes more vertical. However, unlike in the intermediate 
disc case these oscillations are not coupled between the two 
hemispheres. The held lines in the top and bottom hemi¬ 
spheres oscillate and rearrange themselves largely indepen¬ 
dently of one another. For denser discs, poloidal held lines 
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going through the disc are approximately frozen-into the 
disc which moves relatively slowly. Consequently, the field 
motion above and below the disc are uncoupled, allowing 
the field lines to have greater curvature around the mid¬ 
plane of the disc. This is consistent with an observation 
made by Tzeferacos et al (2009), where simulations with 
discs with larger plasma (3 supported field lines with greater 
curvature. These early time asymmetries last on time scales 
of order T ~ 100 ~ 3.5y. In constrast, at late times asym¬ 
metries last on order of T ~ 10 ~ 0.35y. This suggests that 


VzW. ■ 




as the disc has a chance to reach a more stable configura¬ 
tion, timescales of asymmetries will decrease. 


3.2 Stellar Dipole 

A rotating star with a dipole field in the propeller regime 
exhibits quasi-episodic, one sided outflows (Lovelace et al. 
2010). Matter builds up at the magnetospheric boundary, 
compressing the dipole until there is a symmetry breaking 
and the built up matter accretes to the star in a one-sided 
funnel flow. The funnel flow leads to a one-sided outflow in 


Figure 7. Intermediate density dise: The bottom panel shows 
the mass flux M through the upper (solid red line) and lower 
(dashed blue line) boundaries. The middle panel shows the mass 
flux difference, normalized to average total mass flux AM = 
{Mt — Mi))/{Mt + Mb) cis a function of time. The upper plot 
shows the average velocity of outflows, weighted by mass along 
the upper (red) and lower (blue) boundaries. Though the mass 
flux is symmetric during the oscillatory period 100 < t < 400 
the velocity of outflows is greater across the lower boundary. 



Figure 8. Intermediate density dise: Density p (color), poloidal field lines T (white) and coronal mass flux pvp (yellow). The left-hand 
panel is at t = 200, and the center panel is at t = 240. The right-hand panel shows a superposition of the field lines at these two times 
which corresponds to a maxima in the mass outflow in the upper and lower hemispheres, respectively. Note that the field line structure 
exchanges between the two hemispheres at the two times, with the more inclined field lines being associated with stronger outflows. 
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the opposite hemisphere, after which the matter pressure is 
decreased, the magnetosphere re-expands and the process 
can begin anew. 

To minimize the asymmetries induced by disc oscil¬ 
lations, this set of simulations uses the parameters for a 
diffuse disc. However, the simulations with an intermedi¬ 
ate disc showed no qualitative differences in the results. We 
consider the simple case of a stellar dipole aligned with the 
rotation axis of the star. In addition, there may be a disc- 
field. There are two possible configurations - one where the 
two fields are parallel and the other where they are anti¬ 
parallel. 

In the anti-parallel configuration there is an X-type 
magnetic neutral line in the disc (see Fig. 2). As the large 
scale disc-fie Id advects inwards, it annihilates with the 
dipole, allowing matter to accrete onto the star in a con¬ 
tinuous fashion, different from the episodic matter inflows 
of a pure dipole. Despite the additional energy available to 
launch outflows due to the increased stellar accretion rate, 
the held line geometry does not result in enhanced outhows. 

In a parallel conhguration, the large scale field com¬ 
presses the dipole, reducing the magnetospheric radius. The 
key parameters for determining the strengths of magneto- 
spheric outflows is the relative position of the magneto- 
spheric radius Rm and the co-rotation radius Rev When 
Rev < Rm we are in the propeller regime and have strong 
magnet ospheric outflows. Keeping other parameters hxed, 
one can therefore imagine increasing the strength of the 
large disc-held and turning off the propeller. A second 
type of outhow is an inner disc wind, magneto-centrifugally 
launched from open held lines threading the disc. Though 
the disc wind is launched over the entire disc, mass hux 
is dominated by outhows from the inner part of the disc. 


Vz 


-1.0 -0.4 0 0.4 1.0 



Therefore, increasing Rm until it is at the radius of the in¬ 
ner disc wind signihcantly reduces this outhow, as the open 
held lines are pushed outwards. 

For our base case we consider a stellar dipole with di¬ 
mensionless magnetic moment /x = 30. With an anti-parallel 
disc-held of the type considered in the previous section, this 
conhguration initially has an X type neutral line at i? = 6 
(see Fig. [^. The neutral line should be at a radius large 
enough so that the dipole is not fully compressed by the 
disc and its effect essentially negated. However, since most 
outflows from the disc-held are generated at small radii, the 
neutral line should not be too far out either. 

In the parallel held conhguration we consider dipoles 
with jjL — 10, 20, 30. These correspond to initial magne- 
tospheric radii of Rm = 2.8, 4.1, 5 respectively. The co¬ 
rotation radius is Rev = 1.5 so that initially the system is 
in the strong propeller regime. 

3.2.1 Pure Stellar Dipole Field DIP 

This case is useful for comparison with Lovelace et al. (2010) 
and for future runs that include a large scale disc-held. Fig¬ 
ure!^ shows that the outhows are quasi-episodic and dom¬ 
inantly asymmetric. Two types of outhows are observed. In 
the hrst, mass accretes onto the star via a funnel how onto 
either the top or the bottom of the star. At the same time, 
in the opposite hemisphere, the dipole loop expands and re¬ 
connects forming magnetic islands. These magnetic islands 
carry off matter. These were the types of outhows seen by 
Hayashi et al. (1996) for instance, but without the funnel 
how, since their simulations were not bipolar. As in Hayashi 
we see a current sheet in the middle of the expanding loop. 
However, we do not see any signiheant heating of the gas as 


Figure 9. Dense high-jS disc: The bottom panel shows the mass 
hux M through the upper (solid red line) and lower (dashed 
blue line) boundaries. The middle panel shows the mass hux 
difference, normalized to average total mass hux AM = (Mt — 
Mi))/(Mt + Mb) cis a function of time. The upper plot shows the 
average velocity of outhows, weighted by mass along the upper 
(red) and lower (blue) boundaries. 
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Figure 10. Pure Stellar Dipole: Lower panel shows stellar mass 
accretion rate M* (dotted black line line) and mass outflow 
M (red solid line for upper hemisphere; blue dashed line for 
lower hemisphere) as a function of time. The middle panel shows 
the mass flux difference, normalized to average total mass ffux 
AM = {Mt — Mi))/{Mt + Mb) cis a function of time. The up¬ 
per plot shows the average velocity of outffows, weighted by 
mass along the upper (red) and lower (blue) boundaries. The 
right-hand vertical axis for the outffows has been scaled so 
M/M* = 10“^. Note the two types of outffows. When the stellar 
accretion rate correlates to the matter outffow this corresponds 
to funnel ffow accretion in one hemisphere and outffow in the 
other. The three such events correspond to funnel ffow accretion 
in the upper hemisphere and ejections in the lower hemisphere. 
These outffows occur at fast velocity. 

Outffow when there is no stellar accretion corresponds to 
matter pushing against the magnetosphere in one hemisphere 
and being ejected outwards without matter accretion to the 
star. 


the loop expands. These outflows tend to have a greater ve¬ 
locity. In the second type of outflow, the disc pushes on the 
dipole, but has insufficient pressure to accrete onto the star. 
The magnetosphere re-expands, generating an episodic out¬ 
flow (see Fig. ini- Our results are in close agreement with 
Lovelace et al (2010) that found episodic outflows for stellar 
dipoles and Lii et al. (2014) who found episodic outflows for 
stellar dipoles in the strong propeller regime. Both types of 
outflows occur on time scales T ~ 10 ~ 0.35y with a time 
scale T ~ 50 ~ 17y between events. 

Two qualitatively different types of outflows are ob¬ 
served: magnetospheric winds and disc winds. Figure im 
shows a time when both types of outflows occur simula- 
neously. The magnetospheric winds are launched perpen¬ 
dicular to the held lines of the stellar dipole. In contrast, 
the disc wind is launched along open held lines piercing the 
disc and extending off to inhnity. Both outflows are inher¬ 
ently asymmetric, collimated by a toroidal magnetic held 
begining in the disc and extending into the corona. The 
disc wind is driven by held lines that are open in one hemi¬ 
sphere but closed and connected to the star in the other 


Figure 11. Pure Stellar Dipole: Poloidal mass flux p|v| (color), 
magnetic fleld lines (white lines) and coronal poloidal mass flux 
vectors pVp at t = 110 for the case of a pure stellar dipole fleld. 
Note the presence of two outflows, a magnetospheric wind and a 
disc wind. 



Figure 12. Pure Stellar Dipole: The solid red lines are stream¬ 
line of the magnetospheric wind while the solid green lines are 
the streamlines for the disc wind. The dashed green lines are 
magnetic fleld lines anchored at the same point in the disc. All 
data is at t = 110 for a pure stellar dipole fleld with [j = 30. 

hemisphere. This corresponds to an outflow event that was 
not correlated to any stellar accretion. 

Figure shows sample streamlines and held lines of 
the magnetospheric and disc winds along the length of the 
outflow. We note that the magnetospheric wind has a den¬ 
sity 10 times greater than the disc wind, but comparable 
velocity at hxed distance from the launching point. This 
may be useful in explaining observations where jet/counter 
jet densities differ by a factor of a few. Also we note that 
the toroidal magnetic held drops to zero in the center of the 
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Figure 13. Stellar Dipole plus Parallel Disc Field: Poloidal mass 
flux p|v| (color), magnetic fleld lines (white lines) and coronal 
poloidal mass flux vectors p'Vp at t = 87 for the case of a stellar 
dipole and a parallel disc fleld. Note the presence of a magneto- 
spheric and disc wind, though the latter is considerably weaker. 

magnetospheric wind, indicating the presence of a neutral 
layer. 

3.2.2 Stellar Dipole and Parallel Dise-Field PDIP 

These simulations feature two types of episodic outflows. In 
the first, matter builds up at the magnetospheric boundary 
and begins to diffuse across the stellar dipole field. The ini¬ 
tial top-bottom symmetry breaks spontaneously with the 
disc accreting onto the star in a funnel flow in one hemi¬ 
sphere. In the opposite hemisphere, matter is then pushed 
out by the re-expanding magnetosphere owing to the de¬ 
crease in pressure from the disc due to the stellar accretion. 
In this type of ejection the stellar accretion rate and the 
outflow rates correlate strongly. Figure shows a snap¬ 
shot of an episodic magnetospheric and disc wind outflow. 
Figureshows the sporadic nature of these outflows which 
occur in this case at t = 50, 85, 159, & 191. 

The observed outflows are dominated by a magneto¬ 
spheric wind, launched by a re-expanding magnetosphere, 
that had been compressed due to ram pressure from the 
disc. In contrast to the pure dipole the disc wind component 
is suppressed. This is a result of the field line geometry near 
the star. Owing to its centrifugal launching, the disc wind is 
dominated by outflows coming from the innermost region 
of the disc. In the case of a pure stellar dipole field, the 
inner disc supports field lines tilted away from the Z—axis 
which favors the formation of a disc wind. However, in the 
case of a stellar dipole and a parallel disc-field, the field 
lines in the inner disc are nearly vertical. The inclined field 
lines required for a disc wind are pushed out to larger radii 
and the strength of the wind is correspondingly decreased. 
Interestingly, turning down the strength of the disc-field rel¬ 
ative to the dipole field can turn the inner disc wind back 
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Figure 14. Stellar Dipole plus Parallel Disc-Field: Lower panel 
shows stellar mass accretion rate M* (dotted black line line) 
and mass outflow M (red solid line for upper hemisphere; blue 
dashed line for lower hemisphere) as a function of time. The 
middle panel shows the mass flux difference, normalized to av¬ 
erage total mass ffux AM = {Mt — + ^b) cis a func¬ 

tion of time. The upper plot shows the average velocity of out- 
ffows, weighted by mass along the upper (red) and lower (blue) 
boundaries. The right-hand vertical axis for the outffows has been 
scaled so M/M* = 10“^. Note that there are two types of out¬ 
ffows. When the stellar accretion rate correlates with the matter 
outffow this corresponds to funnel ffow accretion in one hemi¬ 
sphere and outffow in the other. These outffows have a greater 
average velocity. Outffow when there is only small accretion cor¬ 
responds to matter pushing against the magnetosphere in one 
hemisphere and being ejected outwards without matter accret¬ 
ing. 

on by allowing for a field line geometry more favorable to 
outflows, as in the case of the pure dipole. 

The parallel disc-field acts to increase the magneto¬ 
spheric wind outflows by increasing the stellar mass flux. 
The loss of angular momentum to outflows acts to increase 
the accretion rate of the disc (Dyda et al. 2012). As a result, 
accretion is increased by roughly a factor of 4, even though 
the midplane magnetic /3 at the magnetospheric boundary 
is roughly contant for the two cases. This increased stellar 
M* results, on average in an increased magnetospheric wind 
mass flux. However, the effect is not great enough to vary 
the time scales of the episodic outflows which in both cases 
remain at roughly order T ~ 50 ~ 17y. This result suggests 
that the addition of a large scale disc-field will increase the 
overall mass flux of outflows, and increase the spread in the 
outflow’s velocity distribution, owing to the fact that the 
magnetosphere is further compressed. 

We can turn off the magnetospheric wind and turn 
back on the inner disc wind by decreasing the strength 
of the dipole field. We consider dimensionless dipoles fi = 
10, 20, 30. At t = 0 these correspond to Rm = 2.8, 4.1, 5, 
respectively. All cases produce strong episodic outflows. 
During such outflows we measure the magnetospheric ra¬ 
dius and find Rm = 1-4, 2.1, 2.2. In the case /x = 10 we are 
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Type 

Rm 

M 


V 

r / Vx {Rm ) 

10 

DIP 

2.6 

0.055 

1.7 

0.48 

0.77 


PDIP 

1.4 

0.009 

0.4 

0.20 

0.40 

20 

DIP 

3.0 

0.078 

0.7 

0.49 

0.85 


PDIP 

1.8 

0.141 

4.2 

0.23 

0.31 

30 

DIP 

3.3 

0.083 

0.1 

0.33 

0.60 


PDIP 

2.1 

0.364 

0.8 

0.38 

0.55 


The largest outbursts were not correlated to mass ac¬ 
creting on the star. Rather there was a re-expant ion of the 
magnetosphere, launching a “knot” of matter. The charac¬ 
teristic density of the “knot” was p ~ 20, roughly 10 times 
greater than magnetospheric ejection for the pure dipole. 
This accounts for the mass flux being a factor of 10 greater. 
We note the velocity of these much stronger outflows are 
roughly a factor of 2 slower than the corresponding mag- 
netespheric ejections in the pure dipole case. 


Table 3. Magnetospheric radius Rrm mass flux M, characteris¬ 
tic velocity v and scaled velocity v/vK{Rm) for different stellar 
dipoles II in the case of a pure stellar dipole DIP and stellar 
dipole and parallel disc-field PDIP. 


no longer in the propeller regime, that is to say Rm ^ Rcr- 
We can compare this to the case of the pure dipole where 
Rm — 2.6, 3.0, 3.3 during outflows. In the case fi — 10, the 
pure dipole outflows are stronger than those with a parallel 
disc-field by a factor of ~ 2 — 3. In contrast, the p = 20, 30 
parallel disc-field cases are stronger than the pure dipole 
cases by a factor of ~ 2 — 5. This points up the importance 
of the field line geometry in determining the strength of 
outflows. Adding magnetic flux in a disc field can hinder 
outflows if it compresses the magnetospheric radius to val¬ 
ues less than Rev On the other hand, if Rm > Rev then 
adding a large scale field can enhance outflows by increas¬ 
ing the accretion rate and providing more energy to drive 
outflows. 

Table gives the magnetospheric radius Rrn, mass out¬ 
flow rate M, characteristic outflow velocity v and scaled 
outflow velocity v/vK{Rm) for different stellar dipole mo¬ 
ments fi in the case of a pure dipole and disc-field during 
typical outflow events correlated with stellar accretion. We 
measure the characteristic outflow velocity by time averag¬ 
ing the velocity of the outflow at Z = 5 during the outflow 
event. All outflows are magnetospheric winds, except for the 
anti-parallel /i = 10 case where magnetospheric ejections 
were absent because Rm < Rev We notice that the mag¬ 
netospheric radius is an increasing function of the dipole 
strength. For constant dipole strength, the addition of a 
disc-field further reduces Rm ■ Further we note that outflow 
mass flux is typically a few percent of the accretion mass 
flux. The outflow velocity is roughly 75% of the Keplerian 
velocity at the launching radius. We note that in the PDIP 
case /i = 10, Rm < Rev- We are no longer in the propellor 
regime and this results in a lower outflow velocity. Increas¬ 
ing /A results in increasing the magnetospheric radius which 
in turn increases the strength of the propellor and the out¬ 
flow velocity. For pure stellar dipoles, the magnetospheric 
radius does not vary appreciably. The dipole is compressed 
until there is an equilibrium between the matter and mag¬ 
netic pressure. This leads to little variation in the size of 
typical outflows. We note that in the case of fi — 10, de¬ 
spite the fact that propeller has been turned off, the pres¬ 
ence of a dipole field does inhibit accretion on the star. This 
induces disc oscillations which leads to asymmetries in the 
outflows. This effect persists until late times, at which point 
the magnetosphere has been completely compressed by the 
advected large scale field and outflows become symmetric. 


3.2.3 Dipole and Anti-Parallel Dise-Field ADIP 


In the case of anti-parallel dipole and disc fields, reconnec¬ 
tion occurs as the disc-field is advected inwards and runs 
into the stellar dipole field. We find funnel flow accretion 
to the poles of the star as found earlier by Miller and Stone 
(1997). This type of accretion is favored for this field geom¬ 
etry because at the X-point there is no magnetic pressure 
and matter can move freely onto the field lines connecting 
the disc and the poles of the star. 

As magnetic flux from the disc-field advects inwards, 
the X-point in the disc moves outwards. This has the effect 
of moving the open field lines threading the disc outwards. 
This causes the characteristic outflow speed, set by the local 
Keplerian velocity, to decrease. Outflows are hence dimin¬ 
ished in this type of field configuration, and should die down 
at late times as sufficient magnetic flux is advected inwards 
and the in ner d isc field lines move outward. 

shows the mass flux M through the upper 


Figure 


15 


and lower hemispheres. As the large scale field is advected 
inwards and annihilates with the dipole, the X-point moves 
into the lower hemisphere. This causes a top/bottom disc 
asymmetry, as matter preferentially falls into the lower 
hemisphere. This leads to a stronger outflow rate in the 
lower hemisphere at late times. Figureshows a snapshot 
of the disc and poloidal magnetic field 

The field line geometry is not favorable for launching 
inner disc winds. Despite initially being in contact with the 
disc (see Fig. [^, field lines anchored to the star quickly 
move into the corona so that there is no efficient mechanism 
for mass loading onto them from the disc. The innermost 
open field lines anchored in the disc are at R > 7, so that the 
disc wind is weak. Based on the velocity cutoff of Vmax ~ 0.2 
in Fig. [T^ we would have estimated this launching radius as 


R?^25. 

The stellar mass accretion rate M* ~ 0.5 which is a 
factor of ~ 10 times smaller than the accretion rate of the 
pure disc-field. In contrast with the pure dipole, annihi¬ 
lations of the disc-field with the stellar dipole reduces the 
dipole strength and allows for matter to more easily accrete 
onto the star. This leads to a nearly constant accretion rate, 
unlike the case of the pure dipole where accretion occurs in 
episodic funnel flow events. The total mass accretion dur¬ 
ing the interval 0 < t < 200 is comparable in the cases 
of the pure dipole and the dipole with parallel disc-field. 
This suggests that the total energy loss in each case is the 
same. The parallel configuration, owing to the lack of out¬ 
flows, converts a much smaller fraction of this energy into 
outflows. 

Figure shows the fraction of mass in outflows with a 
given poloidal velocity for the various combinations of stel- 
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lar dipole and disc field. Over such long time scales the mass 
outflow is to a high degree of accuracy symmetric between 
the two hemispheres so we only plot our results for I’p > 0. 
We further note several interesting features. The distribu¬ 
tions have a sharp velocity cutoff at Vmax = 0.5, 0.65, 1.0 


in the cases of the DIP, PDIP and DISC respectively. These 
velocities correspond to those of a Keplerian disc at radii 
Rk — 1.6, 1.33, 1.0. The velocity cutoff directly corre¬ 
sponds to the Keplerian velocity at the innermost launch¬ 
ing radius. In the case of the disc-field, this is the surface 
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Figure 15. Stellar Dipole plus Anti-Parallel Disc-Field: The 
lower panel shows the mass flux M through the upper (solid 
black line) and lower (dashed red line) boundaries as a function of 
time. The upper panel shows the mass flux difference, normalized 
to average total mass flux AM = {Mt — Mi))/{Mt + M;,) as a 
function of time. We note that outflows are 100 times weaker 
than in the DIP and PDIP cases. 
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Figure 16. Stellar Dipole plus Anti-Parallel Disc Field ADIP: 
Poloidal mass flux p|v| (color), magnetic field lines (white lines) 
and coronal poloidal mass flux vectors pvp at t = 120 for the 
case of a stellar dipole and an anti-parallel disc field. We note 
the weakeness of the disc wind despite the funnel flow accretion 
onto the star. 
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Figure 17. Fraction of total mass outflowed as a function of 
poloidal velocity for DISC (solid black line), DIP (red dashed 
line) , PDIP (blue dash-dot-dot line) and ADIP (green long dash 
line) simulations. The velocity cutoff is set by the Keplerian ve¬ 
locity scale at the innermost launching radius. In the case of 
the disc-field, this is at the surface of the star and therefore 
Vmax = I. In cases with a dipole held this is the magnetospheric 
radius Rm > 1 and therefore r’max < I- For the PDIP case, the 
disc held compresses the magnetosphere and the peak appears 
at larger Vp than for the DIP case. We show the plot for Vp > 0 
since the full distribution is symmetric. 

of the star and therefore Vmax = 1. In cases with a dipole 
field this is the magnetospheric radius Rm > 1 which leads 
to r’max < 1. In the case of the parallel disc-field, the disc- 
field compresses the dipole field and decreases the magne¬ 
tospheric radius. Hence the maximum observed velocity is 
higher than in the case of a pure dipole. Also the peak 
in the distribution is shifted to higher velocities. This pro¬ 
vides information about the field structure near the star: A 
large scale field will have outflows over the entire velocity 
range whereas the presence of a dipole may have a cutoff at 
^max oc Rm^^‘^ siucc outflows will uot bc launched at radii 
R < Rm- 

4 CONCLUSION 

Here we summarize some of our main conclusions. 

1. The distributed disc-fields (DISC) in the absence of a 
stellar field produce symmetric outflows for low density, low 
midplane-/3 discs. This limit agrees with study by Fendt 
and Sheikhnezami (2013). For larger disc densities, there is 
a decoupling of Alfven waves between the upper and lower 
hemispheres and this leads to oscillations in the outflows 
between the two hemispheres. For even higher disc densi¬ 
ties there the jet outflows are strongly asymmetric with the 
outflow changing episodically from one hemisphere to the 
other. 

2. For a stellar dipole field in the absence of a disc-field 
(DIP), two types of outflows are observed - magnetospheric 
winds from the magnetospheric radius Rm and disc winds 
from somewhat larger radii. Both components have similar 


velocities but the magnetospheric winds have significantly 
higher densities and thus carry most of the mass outflow 
M. The asymmetry of the outflows between the top and 
bottom hemispheres increases as the magnetospheric radius 
Rm increases beyond the corotation radius Rcr as found 
earlier by Lovelace et al. (2010). 

3. Introducing a disc-field parallel to the stellar dipole field 
(PDIP) acts to reduce the magnetospheric radius Rm and 
changes the inner disc field structure. The field lines in the 
inner disc field lines are more aligned with the Z—axis and 
as a result the disc wind component is suppressed. At the 
same time the characteristic outflow velocity of the magne¬ 
tospheric winds is increased owing to the reduced magneto¬ 
spheric radius. Outflows are correspondingly more powerful. 
The outflows are episodic and strongly asymmetric about 
the equatorial plane. The asymmetry of the outflows is not 
affected by the initial midplane plasma /3 of the disc. 

4. For a disc-field which is anti-parallel to the stellar dipole 
field (ADIP), we observe reconnection between the field 
components analogous in some respects to the reconnec¬ 
tion between the subsonic solar wind magnetic field and 
the earth’s dipole field at the magnetopause which is inside 
the Earth’s bow shock (e.g., Hargreaves 1992). This recon¬ 
nection acts to prevent a build up of poloidal magnetic flux 
outside the magnetosphere. As a result magnetospheric out¬ 
flows are suppressed. At the same time accretion to the star 
is enhanced. 

Our results suggest that multi-component outflows, 
disc winds and magnetospheric winds, occur for a stellar 
dipole field in the absence of a large-scale disc-field. In¬ 
ward advection of a large scale disc-field acts to turn off 
the disc wind component so that there is only a magneto¬ 
spheric wind. We observe a cutoff in the maximum velocity 
of outflows related to the inner launching radius. This sug¬ 
gests that by measuring the velocity of outflows, one can 
estimate the launching radius. Measurement of the dipole 
component of the magnetic field would then allow one to 
estimate the matter pressure at the magnetospheric bound¬ 
ary. 

Ellerbroek et al. (2013) found they could reproduce the 
knotted structure of outflows from HH1042 with a simple 
periodic outflow model. One possible solution that was con¬ 
sistent with observations was one where the jet and counter 
jet were out of phase by a half period. Observations were 
consistent with a period T ~ 10 — lOOy Our DISC simula¬ 
tion with intermediate jS enters into a regime with period 
T = 2.82y (i^o/O.lAU)^^^. This agrees with the observed 
periods with a reasonable tuning of the inner disc length 
scale Rq. 

Woitas (2002), Hartigan & Hillenbrand (2009) and 
Melnikov et al. (2009) have found that RW Aur A has a jet 
with velocity roughly 2 times greater and density roughly 
2-3 times as large as its counter jet. This is the type of be¬ 
haviour seen in our intermediuate 13 DISC run where in the 
oscillatory phase, the bottom jet mass flux peaks at a rate 
roughly 3 times greater than the top jet and has a velocity 
2-3 times greater. 


(c) 2015 RAS, MNRAS 000,(1^1^ 





16 S. Dyda et al. 


ACKNOWLEDGMENTS 

Resources supporting this work were provided by the NASA 
High-End Computing (HEC) Program through the NASA 
Advanced Supercomputing (NAS) Division at the NASA 
Ames Research Center and the NASA Center for Compu¬ 
tational Sciences (NCCS) at Goddard Space Elight Cen¬ 
ter. The research was supported in part by NASA grant 
NNX14AP30G, and NSE grant AST-1211318. We thank 
the anonymous referee for their very thorough review of 
our work. 


REFERENCES 

Bacciotti, F., Eisloffel, J., & Ray, T.P. 1999, A&A, 350, 917 
Bisnovatyi-Kogan, G.S., & Ruzmaikin, A.A. 1976, Ap&SS, 42, 
401 

Caratti o Garatti, A., Garcia Lopez, R., Weigelt, G., Tambovt- 
seva, L.V., Grinin, V.P., Wheelwright, H., Ilee, J.D., 2013, 
A&A, 554, A66 

Ghagelishvili, G. D., Bodo, G., & Trussoni, E. 1996, A&A, 
306,329 

Goffey, D., Bacciotti, F., Woitas, J., Ray, T.P., & Eisloffel, J. 
2004, ApJ, 604, 758 

Donati, J.-F.;, Jardine, M. M., Gregory, S. G., Petit, P., Bouvier, 
J., Dougados, G., Menard, F., Gollier Gameron, A., Harries, 
T. J., Jeffers, S. V.;, & Paletou, F. 2007, MNRAS, 380, 1297 
Donati, J.-F., Jardine, M. M., Gregory, S. G., Petit, P., Paletou, 

F. ; Bouvier, J., Dougados, G., Menard, F., Gollier Gameron, 
A., Harries, T. J., Hussain, G. A. J., Unruh, Y., Morin, J., 
Marsden, S. G., Manset, N., Aurire, M., Gatala, G., & Ale- 
cian, E. 2008, MNRAS, 386, 1234 

Donati, J.F., Gregory, S.G., Alencar, S.H.P., Hussain, G., Bou¬ 
vier,Dougados, G., J., Jardine, M.M., Menard, F., Romanova, 
M.M.2012, MNRAS, 425, 2948 

Donati, J.F., Gregory, S.G., Alencar, S.H.P., Hussain, G., Bou¬ 
vier, J., Jardine, M.M., Menard, F., Dougados, G., & Ro¬ 
manova, M.M. 2013, MNRAS, 436, 881 
Dyda, S., Lovelace, R.V.E., Ustyugova, G.V., Lii, P.S., , Ro¬ 
manova, M.M., & Koldoba, A.V., 2012, ApJ, 432, 127 
Dyda et al., 2015, in preparation 

Ellerbroek, L. E., Podio, L., Kaper, L., Sana, H., Huppenkothen, 

D., de Koter, A., & Monaco, L. 2013, A&A, 551, A5 
Fendt, G., & Sheikhnezami, S., 2013, ApJ, 774, 12 
Guan, X., & Gammie, G.F. 2009, ApJ, 697, 1901 
Hargreaves, J.K. 1992, The Solar-Terrestrial Enviroment, (Gam- 
bridge Univ. Press: Gambridge), chap. 5 
Hartigan, P., Hillenbrand, L., 2009, ApJ, 705, 1388 
Hayashi, M.R., Shibata, K., & Matsumoto, R. 1996, ApJ, 
468:L37 L40 

Illarionov, A.F., & Sunyaev, R.A., 1975, A&A, 39,185 
Jardine M., Gollier Gameron A., & Donati J.-F. 2002, MNRAS, 
333, 339 

Koldoba, A., Ustyugova, G., Romanova, M.M., Lii, P., Gomins, 

G. , Dyda, S., Lovelace, R.V.E., 2015, New Astronomy, sub¬ 
mitted. 

Lesur, G., Longaretti, P.Y., 2009, A&A, 504, 309. 

Lii, P.S., Romanova, M.M., Ustyugova, G.V., Koldoba, A.V., & 
Lovelace, R.V.E. 2014, MNRAS, 441, 86 
Lii, P., Romanova, M.M., & Lovelace, R.V.E. 2012, MNRAS, 
420, 2020 

Liu, G., Shang, H., Pyo, T., Takami, M., Walter, F., Yan, G., 
Wang, S., Ohashi, N., Hayashi, M., 2012, ApJ, 749, 62. 

Livio M., 1997, in Wickramasinghe D. T., Bicknell G. V., & 
Ferrario L., eds, lAU Goll. 163. ASP Gonf. Ser. Vol. 121; 


Accretion Phenomena and Related Outffows. Astron. Soc. 
Pac., San Francisco, p. 845 

Lovelace, R. V. E., Romanova, M. M., & Bisnovatyi-Kogan, G. 

S. 1999, ApJ, 514, 368 

Lovelace, R. V. E., Romanova, M. M., Ustyugova, G. V., & 
Koldoba, A. V. 2010, MNRAS, 408, 2083 
Lovelace, R.V.E., Romanova, M.M., Lii, P., & Dyda, S. 2014, 
Gomp. Astrophys. and Gosmol., 1, id3, lOpp 
Matsakos, T., Vlahakis, N., Tsinganos, K., Karampelas, K., 
Sauty, G., Gayatte, V., Matt, S. P., Massaglia, S., Trussoni, 

E. , & Mignone, A. 2012, A&A, 545, A53 

Melnikov, S., Eisloffel, J., Bacciotti, F., Woitas, J., Ray, T.P., 
2009, A&A, 506, 763. 

Miller, K.A., Stone, J.M., 1997, ApJ, 489: 890-902. 

Mundt, R., Buehrke, T., Solf, J., Ray, T. P., & Raga, A. G. 1990, 
A&A, 232, 37 

Murphy, G.G., Ferreira, J., Zanni, G., 2010, A&A 512, A82. 
Perrin, M., Graham, J. R., 2007, ApJ, 670, 499 
Ray, T. P., Mundt, R., Dyson, J. E., Falle, S. A. E. G., & Raga, 
A. G. 1996, ApJ, 468, L103 

Ray, T., Dougados, G., Bacciotti, F., Eisl ffel, J., & Ghrysosto- 
mou, A. 2007, Protostars and Planets V, 231 
Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & 

Lovelace, R. V. E. 2004, ApJ, 616, L151 
Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & 

Lovelace, R. V. E. 2005, ApJ, 635, L165 
Romanova, M.M, Ustyugova, G.V., Koldoba, A.V., Ghechetkin, 
V.M., Lovelace, R.V.E. 2002, ApJ, 578, 420 
Sauty, G., Gayatte, V., Lima, J. J. G., Matsakos, T., & 

Tsinganos, K. 2012, ApJ, 759, LI 
Shakura, N.L, & Sunyaev, R.A. 1973, A&A, 24, 337 
Stepanovs, D., Fendt, G., Sheikhnezami, S., 2014, ApJ, 796,29. 
Stute, M., Gracia, J., Vlahakis, N., Tsinganos, K., Mignone, A., 
& Massaglia, S. 2014, MNRAS, 439: 3641 
Tzeferacos, P., Ferrari, A., Mignone, A., Zanni, G., Bodo, G., 
Massaglia, S., 2009, MNRAS, 400, 820. 

Tzeferacos, P., Ferrari, A., Mignone, A., Zanni, G., Bodo, G., 
Massaglia, S., 2013, MNRAS, 428, 3151. 

Ustyugova, G.V., Koldoba, A.V., Romanova, M.M, Lovelace, 
R.V.E. 2006, ApJ, 646, 304 

Velazquez, P.F., Riera, A., Raga, A.G., Toledo-Roy, J.G., 2014, 
ApJ, 794, 128. 

Vidotto, A.A., Gregory, S.G., Jardine, M., Donati, J.F., Petit, 
P., Morin, J., Folsom, G.P., Bouvier, J., Gameron, A.G., Hus¬ 
sain, G., Marsden, S., Waite, LA., Fares, R., Jeffers, S., do 
Nascimento, J.D., 2014, MNRAS, 441:2361-2374. 
von Rekowski, B., Brandenburg, A., Dobler, W. & Shukurov, A. 
2003, A&A, 398, 825 

von Rekowski, B., Brandenburg, A., 2006, Astron. Nachr, 327, 
53. 

Wang, J. G., Sulkanen, M.E, Lovelace, R.V.E., 1992, ApJ, 390, 
46 

Whelan, E.T., Bonito, R., Antoniucci, S., Alcala, J.M., Giannini, 

T. , Nisini, B., Bacciotti, F., Podio, L., Stelzer, B., Gomeron, 

F. , 2014, A&A, 565, A80 

White, M.G., Bicknell, G.V., McGregor, P.J., Salmeron, R., 2014, 
MNRAS, 442, 28. 

Woitas, J., Ray, T.P., Bacciotti, F., Davis, G.J., & Eisloffel, J. 
2002, ApJ, 580, 33 

Zanni, G., Ferrari, A., Rosner, R., Bodo, G., Massaglia, S., 2007, 
A&A 469, 811828 

Zinnecker, H., McGaughrean, M. J., & Rayner, J. T. 1998, Na¬ 
ture, 394, 862 


© 2015 RAS, MNRAS 000, [TUT^ 


